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Abstract 



Soft pomeron exchange at high energy includes minijet production at central rapidity values. This 
predicts copious production of minijets in deep inelastic scattering at small x, which is closely related 
to the fact that vWi exhibits Regge behaviour. It helps to explain also how the integrated inclusive 
cross-section for minijet production in jp or pp collisions can be much greater than the total cross- 
section. 

1 Introduction 

Regge theory successfully predicted M the jp total cross-section measure dP] at HERA. There were 
some other predictions that were much larger^], and it is interesting to discuss why they turned out to 
be wrong. They were based on perturbative calculations of inclusive minijet production. Minijets are 
jets of rather small transverse momentum, say less than 1% of the total available energy y/s. While 
their transverse momentum qx has to be large enough for perturbation theory to be used to predict 
their production rate, it may be too small to allow a clean experimental analysis. Since the scale 
of nonperturbative effects is typically 1 GeV, we would expect that lowest-order perturbation theory 
reproduces the inclusive cross-section da/dqx reasonably successfully down to q™ in ~ 1 GeV, with a 
A'-factor somewhere between \ and 2. One calculates da/dqx from the familiar diagram of figure 1, 
which involves the structure functions of the incoming particles and a hard scattering and yields a pair 
of minijets of approximately equal and opposite transverse momenta. If we integrate it, we obtain M 



where p(q™ in ) is the fraction of events containing minijets, and n is the average number of minijet 
pairs in these events. The naive expectation from figure 1 is that exactly one minijet pair is produced, 
that is n = 1 in (1.1). However, one finds if one evaluates the integral in (1.1) and compares with the 
HERA measurements, that this would require p(q™ in ) > 1, which is wrong by definition. 

So it must be that n > 1. This is not really a surprise. One way in which n can be greater than 1 
is from multiple parton-parton scatterings^] involving two or more partons from each of the initial 
particles. We would expect this to be very important in high-energy nucleus-nucleus reactions^], but 
we doubt whether it is the main mechanism in jp or pp collisions. 

Rather, there is another mechanism W, which is intrinsically nonperturbative. Further, while maybe 
the effects of multiple parton collisons can be approximately handled by an eikonalisation procedure, 
eikonalisation is not relevant for this other mechanism, which may be understood as follows. In 
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Figure 1: Diagram for single-inclusive minijet-pair production 

order to calculate da/dqx from figure f, one needs the two structure functions down to fractional- 
momentum values x of order x ~ q™ in /\/s, that is very small values when ^/s is large. But, from 
elementary kinematics, when a parton of momentum k is pulled out of a particle of momentum p with 
k = xp + . . . , the squared invariant mass of the system it leaves behind is^] 

so 1.2 

X 

and so is large when x is small. That is, the upper and lower clusters of residual fragments of the initial 
hadrons in figure f each have large invariant mass and so they are very likely to contain additional 
minijets in some fraction of the events. 

Of course, the structure function is measured directly in deep inelastic lepton scattering, where also 
the residual fragments of the proton will have large invariant mass at small x and so are likely to 
include minijets. These minijets will be distributed uniformly in rapidity in the central region. In 
this paper we calculate this effect. We expect it to be of significant magnitude, because it is closely 
connected with the reason for structure functions exhibiting Regge behaviour at small x. 




Figure 2: Parton model for vWi 



Consider vWi at a value of Q 2 that is only moderately large, so that perturbative evolution has not 
yet set in and the parton model applies: figure 2. The squared energy of the lower bubble, which is the 
amplitude for finding a parton k in the proton p, again satisfies (1.2). Since this bubble is an elastic 
strong-interaction amplitude, its high-energy behaviour is governed by Regge theory and is just a sum 
of terms Sq' '. If we insert this into the calculation of vWi from figure 2, we obtain^] a sum of terms 
(l/x)"^- 1 . NMC data[ § ] at moderate Q 2 and not-too-small x show that vWi contains such Regge 
terms at small x, approximately a constant from soft-pomeron exchange and close to ^fx from fn a i 
exchange^]. 
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it is not widely Known tnat sucn a Denaviour, wnicn we stress is nonperturoative, arises 1 - J Decause 
the proton fragments that remain when a small- a; parton is pulled out have large invariant mass. 
We note that the Monte Carlos used at HERA set both k 2 and to zero (before the perturbative 
evolution begins) and so do not make so large. Thus they miss a nonperturbative effect which is surely 
important. 

The occasional production of high-py jets is part of the Regge behaviour; in particular, it is part 
of soft pomeron exchange. As the energy increases, this minijet production becomes more and more 
frequent, but it is still part of soft pomeron exchange and does not change the effective power of s with 
which the total cross section rises. We calculate the probability that soft pomeron exchange includes 
such production. 

In section 2, we introduce the coupling of the soft pomeron to a gluon and estimate its strength from 
the small- a; behaviour of the proton's gluon distribution. We use this in section 3 to estimate minijet 
production in deep inelastic lepton scattering at small x, and find that it should be copious for small 

~min 
q T . 

2 Inclusive cross-section for production of a single pair of minijets 

As is well known^], at small x there is no k^ ordering, so that perturbative and nonperturbative 
effects are inevitably intermingled. In particular, as the dominant mechanism for the production of 
two pairs of minijets in hadron-hadron collisions, we propose figure 3. As we go down the diagram, 
soft pomerons alternate with hard interactions that produce the minijets. If we slice the diagram down 
the middle, the left half corresponds to figure 1, with the upper (or lower) hard interaction in figure 

3 corresponding to the one shown explicitly in figure 1, and the other one being implicitly included 
within the lower (or upper) bubble in figure 1. Gluonic partons are more efficient at producing high-py 
jets than are quarks^], so we concentrate on the case where the hard collisions are glue-glue. 

This means that we need the coupling of the soft pomeron to a gluon. In order to extract this, we first 
calculate the inclusive cross-section for production of a single pair of minijets. That is, we calculate 
figure 4. For the coupling of the pomeron to an on-shell quark we have used previously^ 

M a (2-la) 

together with a signature factor that ensures that quarks and antiquarks couple with equal sign. This 
is very successful phenomenologically, with 

/3 ~ 2 GeV" 1 (2.2a) 

We need the coupling of the soft pomeron to a gluon only at zero momentum transfer t. At t = the 
most general coupling to a gluon of momentum k would be 

S ab A(k 2 ) <f v k a (2.16) 

together with terms in k^ or k v (or both), where fj, and v are Lorentz indices for the gluon and a and 
6 are the corresponding colour indices. To construct a high-energy forward scattering amplitude for 
either qq, qg or gg scattering, we multiply a corresponding pair of couplings (2.1) by (a'so) e , where 
so is the squared centre-of-mass energy for the amplitude and e and a' are the intercept and slope of 
the pomeron trajectory: 

e « 0.08 a' = 0.25 GeV" 2 (2.26) 
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Figure 3: Diagram for double-inclusive minijet-pair production, with alternating soft pomeron exchange and 
hard scattering 



When the gluon is on-shell, only the term in its coupling shown explicitly in (2.1b) contributes, as the 
others vanish when multiplied by the gluon polarisation vectors. In most applications where the gluon 
is not on shell, it is still true that we need only this term, because the other ends of the gluon lines 
are coupled to an amplitude that is annihilated by k^ or k v . We include in A(k 2 ) a (nonperturbative) 
propagator for each gluon leg. As we shall see, we do not need to know the function A(k 2 ), but only 
the value of the integral 

9 = ^ (2 + e)- 1 J ^ dk 2 (-k 2 ) 2 +*A(k 2 ) (2.3) 

The calculation of figure 4 involves the square M of the hard-scattering gg amplitude, integrated 
over the momenta of the final-state jets - or rather, because of momentum conservation, over the 
momentum q of just one of them. Because of the g 111 '5 ab in the coupling (2.1b) of the pomeron to a 
gluon, the Lorentz indices and colour indices for M are summed in the same way as in the calculation 
of the differential cross-section da/dt for the sub-process. In fact, 

512s 2 f d 3 q ^da . 

M = / -±6(s + t + u)-p: 2.4 

7T J q° dt 

To calculate figure 4 we introduce momenta k\ and ki as indicated, and express each in the form 

ki = Xip - mp' + k iT (2.5) 

In our calculation, we may neglect the masses of the particles p and p' . Although there is no kx 
ordering, there is x and y ordering 

[10] 

in the region of phase space that dominates in the calculation: 
1 > x x > x 2 > 0, < yt < t/2 < 1 (2.6) 
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Figure 4: Diagram for single-inclusive minijet-pair production, with hard scattering between soft pomeron 
exchanges 



What we mean by these strong inequalities is that if we write x\ ~ (a's) Vl and x 2 ~ (a's) V2 then 
< i]i < i]2. We use 

d 4 ki ~ Itt / — dk\ dk 2 1T [ d 4 k 2 ~ lyr / ^rffc 2 2 rfk^ T (2.7a) 



"1 ~ 7 J Vl 

The subenergies across the two pomerons satisfy 

A/i "I - ki A/o ~~ I - ^2T 1 / 7\ 

s i ~ Ui s ~ s 2 ~ X 2 S ~ (2.7o) 

•''I ~ V2 

The subenergy si is raised to the power e and is multiplied by the coupling of the gluon to the pomeron, 
so that we encounter the integral 

, -*? 

2 / .71,2 / ,2 ,2 + 6 



(2tt) 



difcf / dki T (-kf - \t{ T y +e A{k[) (2.7c) 



where the upper limit on the k\x integration arises from the need to keep si positive. This is the only 
way in which the gluon coupling enters the calculation; on performing the k\x integration we retrieve 
(2-3). ' " 

The coupling function A(k 2 ) of the gluon to the soft pomeron is nonperturbative, and we have defined 
it to include also propagators for the two gluon legs. Thus we expect 

[12] 

that it is appreciable in size 

only when \k 2 \ is less than about 1 GeV 2 . We see from (2.7c) that there is a similar constraint on ky, 
so the partons k\ and k 2 have only rather small transverse momentum and the transverse momenta of 
the two minijets produced in the hard collison are almost equal and opposite. So, knowing that each 
jet is paired with a second one in this way, we calculate the inclusive cross cross section for producing 
a pair of minjets: 

da pail 1 f , , , ,si2 da 



- - / dx 1 dy 2 F(x 1 )F(y 2 )S(s 12 + t 12 + u l2 ) — — (2.8a) 

a°q z J 7r at\ 2 
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wiin 6i2,t-i2 anu "12 Lne iviancieisiam vanaoies 101 me cenirai narci scauering ^so inai s'12 = {K\ — 
k 2 ) 2 ~ xiy 2 s), and 

F(x) = 24tt/3 gx~ l - e (2.9) 
From the familiar structure of (2.8a), we identify F(x) as the gluon structure function of the proton 
at small x. From its known normalisation^], we deduce that for Q 2 ~ 5 GeV 2 , 



g « 15 MeV 

If we integrate over all qx > q™ in , M in (2.4) becomes 



M(s) = 512s / df-^0(s 



4q 



min 

T 



and the integrated inclusive cross-section for producing a pair of minijets is 



dqi 



97r 2 /3 y 

16s 



dx 1 dy 2 (x 1 y 2 ) 



-2-E 



Ml 



(2.10) 



(2.11) 



(2.12) 



We change integration variables from y 2 to s\ 2 
to (2.6), which has become 



:\y 2 s. We then perform the x\ integration, subject 



Si 2 

1 > zi > — 



We also obtain from (2.6) the conditions S12 ^> —(k 2 + k 2 



(2.13) 

? T ) (i = 1,2), but these can be ignored 
because (2.11) requires S12 > 4g™ in , and the form factor that couples the soft pomeron to the gluon 
suppresses contributions to the integral (2.3) from large (k 2 + k 2 T ). Because the limits (2.13) on the 
x\ integration only enter in a logarithm, at high energy it is a good approximation to replace the 
strong conditions (2.13) by simple inequalities. Thus, with o TOT = 18/3o(a's) e , (2.12) becomes 



h 



rrTOT , 



dqi 



da pa 



dqj 



32a' e 



rfsi 2 s 12 2 e M(si 2 )log 



4g5? ln 



S\2 



(2.14) 



where M{s\ 2 ) is given in (2.11). 

Although we have explicitly considered pp scattering, the result (2.14) is valid for the scattering of 
any pair of hadrons. It is valid also when one (or both) of the hadrons is replaced with a real photon. 
This is because, in minijet production, as in the total cross section, the photon behaves like a hadron: 
the "direct" component of its structure function is unimportant^]. 
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Figure 5: The integral I\ of (2.14) 
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vve are unaoie 10 caicuiaie y^-i-i) anaiy ucany. iNumericai compuiauon, uaKing accouni oniy 01 gg — f gg 
in the hard scattering, gives the results shown in figure 5. We have used a running coupling a s 
evaluated at 4g™ in and with Aqcd = 250 MeV, and introduced a Ji'-factor of 2 to take account of 
non-leading-order effects. We recall that the integrated inclusive cross section has the interpretation 
given in (1.1), so that it is no surprise that it is larger than the total cross-section. 

As a check on our calculation, we have compared our results with the UA1 minijet data^] in figure 
6. In this experiment the jet event cross section was defined as the cross section for producing at 
least one cluster with E™ w >5 GeV at any rapidity. A nominal 5 GeV cut on the cluster transverse 
energy is quoted as corresponding to an effective average parton px threshold of between 3 and 4 GeV 
and this effective threshold depends, in an unspecified way, upon the cm. energy. We have chosen 
to compare our results with the data using a px cut of 4 GeV, and as can be seen from figure 6 the 
comparison is reasonable. As our calculation is asymptotic, in the sense that we have taken the x — > 
limit (our structure function is simply x~ e , with no accompanying power of (1 — x) ), our estimate of 
the inclusive minijet cross section is less meaningful at the lower energies than at the higher. 




200 



900 



Figure 6: Comparison of calculated inclusive minijet cross-section integrated down to qx 
datat 15 ] 



4 GeV with UA1 



3. Deep inelastic lepton scattering at small x 

We now calculate how much minijet production may be expected in the central rapidity region in deep 
inelastic scattering. 

Pomeron exchange is supposed to dominated in the small- a; behaviour of vWi- At small and moderate 
values of Q 2 it is expected that the pomeron that is involved is the soft pomeron, so that 

vW 2 ~ C{Q 2 )x~ e (3.1) 

with e ~ 0.08 as is given in (2.2b). This expectation is well verifie d^ by NMC data[ § ] at moderate 
Q 2 and not-too-small x, though the measurement at Q 2 = 8.5 GeV 2 by the HI collaboration^] gives 
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a larger enecuve vame 01 t ai very smaii x ana raises me liiieresiing possioiniy inai someining eise is 
involved, maybe even the nonperturbative pomeron 

[10] 

Even if the effective value of e increases with Q 2 , the effect of this should to some extent cancel in the 
ratio of the minijet cross-section to the total cross-section (see (3.3) below). In any case, whatever 
mechanism is responsible for the rise in the effective value of e will surely increase minijet production. 
So we here confine our attention to the contribution from soft pomeron exchange, and interpret our 
results as lower limits for the fraction of events that contain minijets. We assume that soft pomeron 
exchange factorises, so that (3.1) corresponds to figure 7a. The probability that the soft pomeron 
exchange includes the production of a pair of minijets is then calculated from figure 7b. As before, 
we expect that the main mechanism for minijet production is gluon-initiated, and confine ourselves to 
considering this. 




(a) (b) 



Figure 7: (a) Parton model for vW 2 at small a;16(6) Single-inclusive minijet-pair production in deep inelastic 
scattering at small x 

The calculation is similar to that outlined in section 2, with p' replaced with Q and the condition on 
X2 in (2.6) strengthened so that x 2 x. This has the consequence that (2.13) becomes 

XSi 2 

1 > x l > — ^ 3.2 
I 1 

where /2 2 = — {k\ + k?, T ). Strictly, the appearance of /2 2 here prevents us from factoring off the k\ 
and k?, T integrations and absorbing them into the integral (2.7c) that defines g, but at very large s 
the error in this is small. This is because, as we have explained, the pomeron's coupling to the gluon, 
together with the nonperturbative gluon propagators, effectively constrain ^ to be of the order of 1 
GeV or so, and since it will appear only in a logarithm we can set it constant without too much error. 
We find that, for Q 2 > /2 2 , the fraction of events in which there is a pair of minijets with qx > q™ in is 

1 r d Tr 2 n 2 f j2 l x f ii 2 \ 

w , n2 , / dq T —vW 2 = -X ds 12 s- 2 - e M(s 12 ) log M^- (3.3) 

One might have thought that the argument of the logarithm here should be Q 2 /xsi 2 , that is s/s\ 2 , 
in conformity with (2.14). At extremely large values of s (or extremely small values of x) the two are 
equivalent, but at subasymptotic values the kinematics lead to the form in (3.3). 

We plot the fraction (3.3) in figure 8 for two values of q™ in . The bands in the figure correspond to 
fi 2 varying in the range 0.5 to 2 GeV 2 . Note that q™ in refers to the parton: contributions from the 
underlying event background must be subtracted from data. 
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Figure 8: Fraction of deep-inelastic-scattering events that contain minijets 
4 Production of two pairs of minijets in pp or jp collisions 

The calculation of figure 3 is along the same lines. It yields the inclusive cross section for the production 
of two pairs of minijets. The transverse momenta of the minijets in each pair are almost equal and 
opposite, and the total-momentum vectors of the two pairs are randomly distributed in rapidity in 
the central region of rapidity, and therefore, since we assume that the total energy is very high and a 
large rapidity interval is available, they are usually well separated. We express each of the four ki in 
the form (2.5), and obtain the dominant contribution from x and y ordering: 

1 > x x > x 2 > x 3 > x 4 > < y x < y 2 < V3 < J/4 < 1 (4.1) 

The integrated inclusive cross-section for producing two pairs of minijets is 

dq T dq' T - — — = n _ ° / dxidy 2 dx ?) dy i {xiy2X ? ,y i s)~ 2 ~ t M{xiy 2 s)M{x ?) y i s) (4.2) 
aqTdq T 2bba' e J 

We change integration variables from y 2 to S12 = x\y 2 s and from £3 to S34 = x^y^s. We then perform 
the x\ and y± integrations, subject to the conditions (4.1). The important conditions read 



?/2>Sl2/s £3>S34/s ?/2«3S<M (4.3a 



where now 

With these conditions, 
with 



= min (-k z 2 - k z 2T , -ki - k z 3T ) (4.36) 

f ^l.ilog^ 9 {R-l) (4.4a) 
Xi J x 3 

R = (4.46) 

•S12S34 

As we have explained, the appearance of fi 2 prevents us from factoring off the , A;§ , k| T and k| T 
integrations and absorbing them into the integral (2.7c) that defines g, At extremely large s the error 
in doing this is small. The result of setting /2 2 constant is that (4.2) becomes 

1 f , d 2 a pa[TS 

h = TOT / dq T dq' T — — — - 
a Jq™ in dq^dq T 
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Figure 9: The integral I 2 of (4.5) 



A, A 



= 1Q 24^, 2e / mm2 ds i2ds 34 (s 12 s 3i )- 2 - e M(s 12 )M(s 3i ) log 2 R 0(R - 1) (4.5) 

In figure 9 we plot this for various values of q™ in . The bands again correspond to p 2 in the range 0.5 
to 2 GeV 2 . 

The interpretation of the doubly-integrated inclusive cross section 3.2) is analogous to (1.1): 

h < n(n - l)p(qT n ) (4.6) 

where again p(q™ in ) is the fraction of events that contain at least one minijet pair, and n is the average 
number of pairs of minijets (including the pair containing the "trigger" minijet) in those events. We 
have written (4.6) as an inequality, in order to account of the other mechanisms for providing an 
inclusive cross-section for two minijet pairs, for example those involving more than one parton from 
each initial particle^]. The need for the presence of some other mechanism is evident because, from 
its definition, p(q™ in ) < 1, so that it must be that 

irh * 1 (4 - 7) 



5 Discussion 

In this paper we have considered inclusive cross-sections for minijet production, and have emphasised 
that they are not related in a simple way to total cross sections. We have found that there is a 
nonperturbative mechanism, which forms part of soft pomeron exchange and is closely related to the 
small- a; behaviour of structure functions, that for small q™ in provides copious minijet production in 
high-energy hadron-hadron, photon-hadron and photon-photon collisions, and in deep inelastic lepton 
scattering at small x. The mechanism is not included in the standard Monte Carlos, for example those 
used to analyse HERA data, and we urge that this be put right. 

The mechanism is equally applicable to the production of heavy flavour in the central region, and 
yields charm and beauty cross sections in high energy hadron-hadron interactions that are compatible 
with data. As is the case for the minijets, this does not affect the overall rise in the total cross section. 
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±ne nonpenuiDauve mecnanism 01 course memoes penuiDauve nara-scauermg suo-processes. vve 
have included only glue-glue subprocesses, but for minijet production the consequent under-estimation 
will be very small. We have also worked only to lowest order in the perturbative as, evaluated at 
4g™ in with Aqcd = 250 MeV. We have rather arbitrarily chosen a Ji'-factor of 2, and there is also 
the uncertainty on the value of g, mainly arising from the poorly-known gluon structure function at 
very small x. Nonetheless for the inclusive minijet cross section the order of magnitude is sensible, as 
is shown by the comparison in figure 6 with UA1 data. 

The calculations show that it is incorrect to ascribe the rise in the photon-hadron or hadron-hadron 
total cross section to the onset of minijet production. The minijets can be seen to be an integral part 
of the pomeron, which is itself the controlling mechanism for the total cross section. 

This research is supported in part by the EC Programme "Human Capital and Mobility" Network 
"Physics at High Energy Colliders" contract CHRX-CT93-0537 (DG 12 COMA) 
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